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We study charge and spin transport through an interacting quantum wire, caused by backscat- 
tering off an effective impurity potential with a periodic time-dependence. The adiabatic regime 
of this pump for charge and spin is shown to depend on the presence of interactions in the wire. 
For the non-interacting case the transported quantities in a period are calculated and found to be 
adiabatic (independent of pumping frequency) for all frequencies £1 <C vf/£, where £ is the range 
of the scattering potential, and vf is the Fermi velocity; this result is along the lines of adiabatic 
pumping in mesoscopic systems. On the other hand, we show that for a wire with repulsive electron- 
electron interactions the adiabatic regime is confined to TiQ. <C cor, with ur an energy scale set by the 
backscattering potential. Using symmetry and scaling properties of the quantum wire Hamiltonian 
we relate the charge Q and spin S transported through the wire in a period 2n/Q. to an integral 
involving quasi-static backscattering conductances. We also show that the pumped charge Q (or 
the spin S) is quantized in the adiabatic limit if the conductance of the system is zero at the stable 
fixed point of the renormalization-group (RG) transformations. The quanta transported are given 
by the winding number of a complex coupling F(fif) as it traverses a closed path C enclosing the 
origin r — 0. No adiabatic charge or spin is transported when the RG fixed point corresponds to 
a perfectly conducting wire, or if the path C does not enclose the point T — 0. By contrast, for a 
RG marginal conductance - which is the case for non-interacting electrons - the charge transported 
in a cycle is not quantized, rather it is proportional to the area enclosed by the path or circuit C, 
in agreement with previous studies of non-interacting mesoscopic systems. Finite size, tempera- 
ture, and frequency effects on the transported quantities follow from the relation between adiabatic 
transport coefficients (backscattering conductances) and the quasi-static conductance of the system. 

PACS numbers: 71.10.Pm, 72.25.Pn, 73.63.Nm, 73.63.Fg 



I. INTRODUCTION 

Quantum adiabatic charge transport through a system 
is achieved by a slow periodic modulation in time of a 
set of external potentials. This phenomenon has been 
investigated in one-dimensional non-interacting systems 
with gapped excitations (band insulators)aau, for which 
case the pumped charge Q (charge transported through 
the system in a single period) was shown to be quantized. 
In the context of adiabatic transport through- quantum 
dots, recent theoretical work on pumpingptrO finds no 
such universal behavior for Q, which strongly depends 
on how the pump is operated. On the other hand, for 
gapless systems in one space dimension (ID Luttingcr 
liquids), it is possible to pump weakly quantized charge 
and spin by the temporal variation of a spatially localized 
potentialQ. 

In our earlier work, Ref. [?|, the weak quantization 
of adiabatically pumped charge or spin in a Luttingcr 
liquid was argued to follow from the repulsive nature 
of electron-electron interactions. These arguments were 
given by solving exactly the special case of the spinlcss 
Luttingcr model with interaction parameter K = 1/2 
in the presence of an impurity with arbitrary time- 
dependence r(J7t). Quantization of pumped charge was 
shown in the limit £1 — > with power law corrections 
AQ ~ (fi/wr) determined by the barrier energy scale 
wr- The quanta transported are given by the wind- 
ing number of the complex backscattering amplitude T. 
Therefore, the pumped quantities arc invariant under 



(small) continuous deformations of the path C, around 
the origin of the complex plane, traced by T. 

In this paper, we prove the adiabatic quantization of 
the pumped charge and spin for gapless interacting sys- 
tems using the properties of chiral symmetry and scale 
invariancc. We derive a connection between the adia- 
batically pumped quantities and an integral involving 
the quasistatic (equilibrium) value of the conductance 
through the system. The quasistatic conductance is cal- 
culated using an instantaneous renormalized (because of 
bulk interactions) backscattering amplitude T(t). For ex- 
ample, in its simplest form for a spinless system, the ex- 
pression for the pumped charge is: 

Q=^--i d ^[^-^Gr\, (1) 
2m J I 

where G r is the equilibrium conductance of the one- 
dimensional wire with the instantaneous (renormalized) 
backscattering amplitude T(t). Our derivation generi- 
cally links quantization of adiabatic charge transport to 
insulating rcnormalization group (RG) fixed points of 
the interacting wire with an impurity. Thus, whenever 
the backscattering potential leads to insulating behav- 
ior (G r — > 0) in response to an applied DC voltage, the 
pumped charge is perfectly quantized, albeit the time to 
transport such charge adiabatically is large (we discuss 
in this paper the frequency scales for adiabaticity) . The 
quanta is given by the winding number (about the ori- 
gin r = 0) of the backscattering amplitude T. For the 
particular case of non-interacting electrons, for which the 
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backscattering amplitude is marginal in the RG sense, the 
aforementioned quantization is lost, and we show that 
our result recovers the "area law" relation obtained by 
Brouwercl. 

We emphasize that the relation between the adiabatic 
charge and spin transported in a periodic cycle and the 
quasistatic conductance is rather general; it applies to 
systems with both short and long-range bulk interactions, 
and is derived here without taking recourse to bosoniza- 
tion. Within the bosonization approach, arguments have 
been put forth in Ref. || that the quantization is a conse- 
quence of pinning the bosonic field at the impurity to a 
minimum of a single cosine potential when the backscat- 
tering potential is relevant; the displacement of the bot- 
tom of the potential after a complete cycle leads to the 
quantization condition. Here we show that the mecha- 
nism of pumping is based on controlling the spatially lo- 
calized backscattering process that breaks a continuous 
global (chiral) symmetry of the bulk one-dimensional sys- 
tem, and point out the crucial role of interactions in the 
wire in determining the transported quantities. 

The rest of this paper is organized as follows. In Sec- 
tion II we consider the case of charge pumping through a 
quantum wire of interacting spinless fermions with gap- 
less excitations. Wc introduce the model in Section IIA 
and relate it to a time dependent single impurity model, 
with a (complex) coupling T(t). In Section IIB we show 
how the temporal variation of the phase of T can lead 
to the transport of charge Q through an otherwise clean 
quantum wire. This is done by making use of the Ward 
identities which arc derived in that section. In Section 
IIC we obtain a relation between the pumped charge in a 
cycle and the backscattering-current conductance of the 
wire with renormalized parameter. In Section IID we re- 
late the pumped quanta to the usual conductance of the 
quantum wire. In Section III we consider electrons with 
spin and derive expressions for the spin (S) and charge 
(Q) transported per cycle by the spin and charge pumps. 
These expressions allow us to look at the limits of quan- 
tization and at the effect of finite temperature and finite 
size on the value of Q and S, which is the content of Sec- 
tion IV. In Section V we solve the time-dependent im- 
purity problem exactly for non-interacting electrons, and 
also present details of the calculations used in our earlier 
publication, Ref. 0, for the spinless Luttinger model at 
the K = 1/2 point. We identify the adiabatic regime, and 
compare the results with the approach of Sections II and 
III. Finally, in Section VI we conclude with a summary 
of our main contributions. 



II. TRANSPORT OF SPINLESS ELECTRONS 
BY PARAMETRIC PUMPING 

A. Model Hamiltonian 

An adequate low energy description of a system of ID 
fermions is obtained by separating the Fermi field ip( x ) 



into right (R) and left (L) moving chiral fields near the 
Fermi points ±kp 



V>(z, t) = e ik ^ip R {x, t) + e- ikp ^ L {x, t). 



(2) 



This turns a non-relativistic action into a relativistic 
Dirac fermion action in (1+1) space-time dimensions, 
withpthc Fermi velocity vf playing the role of the speed of 
lightll Our notation assumes spinless fermions although, 
as we show later, it can be straightforwardly generalized 
to include spin. We consider such a system with a Haejiil- 
tonian H that preserves the global chiral symmetrycJ 



*(a0, 



V(x)e 



2Q7 



where 



IpR 



and * = l-V = (tpt,ipi). 



(3) 



(4) 



Here the 7-matrices are 7 ' 1 = (J x ,y, and 7 = i7°7 1 = 
—<j z , with cr x ,y,z being the Pauli matrices; and we choose 
vp = e = h = 1. The chiral symmetry ensures that 
the axial charge J dx j = J dx is conserved in 

the absence of an external electromagnetic field ("chiral 
anomaly" ) . 

Adding a short range impurity potential V(x) - non- 
zero only in a finite range of length !<L, the length of 
the wire - gives rise to an additional term in the Hamil- 
tonian H: 



dx V(x)tp^ (x)ijj(x) 

dx V(x)e i2kFX ip[(x)^ R {x) + H.c. (5) 
+ j dxV{x) ip R {x)ip R (x) + R -> L . 



Here the first term describes the backscattering of right 
movers into left movers and vice versa, while the sec- 
ond term describes forward scattering of both right and 
left movers with identical scattering amplitude and phase 
shifts. The idea of parametric pumping in a quantum 
wireO is to generate a current by varying parameters that 
control the "shape" of this scattering potential. Since 
the current (axial charge j) involves the difference in the 
number of right and left movers in the wire, the forward 
scattering term, which docs not distinguish between the 
right and left movers, plays no role in generating this cur- 
rent. It is only the backscattering terms which can lead 
to a non-zero current when the potential V(x) is suit- 
ably manipulated. Thus, for the purpose of determining 
the pumped current at low energies (when the continuum 
field theory description (Q) holds), we can write down an 
effective Hamiltonian which describes the most relevant 
(in the RG sense) backscattering processes. Thus the im- 
purity contribution to H can be written as a local (x = 0) 
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FIG. 1: Variation of V due to V(x,t) changing shape by 
changing its value about Ef- 



backscattering term: 
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= TQ+Q^T* 
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(G) 
(7) 



Here Q represents the local operator ipL^niO) and T = 
J dx V(x)e i2kFX represents the 2kp - Fourier-transform of 
the impurity potential. An important point about this 
mapping from a local potential V(x) on to a delta func- 
tion potential T, is that it is valid only at energies lower 
than ^njp/l where £ is a measure of the range of the 
potentials. Let us note here that by simple dimension 
counting (in natural units, with h = e = vp = 1) the 
coupling r has zero dimension. 

If we vary the "shape" of the potential V(x) in time, 
then it amounts to varying both the phase <p and ampli- 
tude of r in time. This can be understood qualitatively as 
follows (see also Fig |l|). Consider a time [t = in Fig [l]) 
when the potential profile is symmetric about x = 0, 
then the Fourier transform is real and the phase </> = 0. 
At a later time, when the profile is anti-symmetric, the 
phase is (f> = n/2, and then, as the potential changes sign, 
<t> = 7r until the potential regains its original shape and 
the phase goes to <f> = 2ir. Thus, by changing the po- 
tential profile, one gets an effective time dependent delta 
function potential |T(t)|e^W 5(x). 



1. Double barriers: a realization ofT(t) 

In the pumping device proposed in Rcf. |^, the changing 
V(x,t) of Fig. [y is realized by oscillating two potential 
barriers placed some distance apart: 



V(x,t) = V -(x)f-(t) + V + (x)f+(t), 



(8) 



with out of phase periodic functions / + and /_ which 
change sign in a cycle. For example, f-(t) = cos cot and 



f + (t) = cos(tot + 9). We identify 



r(t) = V -(2k F )f-(t) + V+(2k F )f + (t), where 



dx e 



V ~' + (x). 



(9) 



If V ' + (x) is asymmetrical about x = 0, then the Fourier 

transforms Vy _ ' + are n °t rea ^ anc ^ as a resu lt I\t) is com- 
plex. Furthermore, if the functions are different so that 
Arg[V^] ^ ArgfVj^] then by varying the amplitudes and 

the signs of the complex quantities f_r + \V (accom- 
plished by varying f±) one can reach any point in the 
complex r plane. Thus, by choosing appropriate func- 
tional forms f± (t) , one can trace any desired path in the 
complex T plane (as in Fig. pi). 



2. Special case: simple form for V(x,t) 

We choose delta function forms for V^~(x) and Vq~(x): 
V Q -(x) = 7 -<5(.t + t/2), V+(x) = j+S{x - 1/2) 
to obtain: 

V(x, t) = j~5(x + 1/2) cos tot + j + 5(x - 1/2) cos(wi + 6). 

This gives an effective coupling 

T = ~f- e - lkFi cos Lot + ~f+e ,kFe cos(Lut + 9), 

which is of the form T(t) = T +U e %ult + T^e" 1 "*, where 
T± w = V + e-' lkpe + V Q ~e ikFe±ie . Written in this way, 
the time dependence of T can be distinguished from that 
obtained by applying a voltage - when only one of T± u 
is present, as in Ref. [l(]. 

The phase 4> of the complex valued parameter V is 
time dependent whenever 2kpt ^ nir and 9 =/= mx 
(n = 0, ±1 . . .). The first condition is the off-resonance 
condition, essential for backscattering to occur, while the 
second is essential for traversing a closed path in the T 
plane with a non-zero area. The rate of change of this 
phase, 4>, can be written as: 



lo^i 7 + sa\9sa\2kpl 



where 



X 2 (t) = cos 2 k F t [y + cos(ojt + 9) +7 cosLot] 2 
+ sin 2 kpl [7 + cos(uot + 9) — 7 ~ coswt] . 

It is clear from this expression that 4>(t) does not change 
sign with time, and thereby takes the complex quantity 
r around a contour, as shown in Fig [I]. 
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B. Pumping parameter for the charge pump 

In the rest of this section we will show that the phase 
of the complex quantity T is the pumping parameter. 
The total Hamiltonian is: 



H = H + H T 



(10) 



The backscattering term (0) breaks the continuous chiral 
symmetry (|J), and thus gives rise to the backscattering 
current 



J r (t) 



1 



dx dtj(t, x) 



(11) 



= - ldx\j(t,x),H T }=iTQ(t)+-a.c, (12) 



which measures the rate of change of total charge of the 
right movers in the system due to the backward scattering 
of the chiral fermions. It is important to note that we do 
not consider, in Ti, the forward scattering term arising 
from the impurity potential. This is because forward 
scattering involves local fermion density and therefore 
conserves the chiral charge. 

In order to look at the equilibrium expectation value 
of (111]) we consider the functional 



Z\r\ = {S T (oo,-ao)), 



(13) 



with a constant (other then the adiabatic switching on 
and off) source term Hr that connects the vacuum state 
of H at times t = — oo with the same state at t = oo. 
As a consequence of the global chiral symmetry a change 
of variables, as in (|||), with an infinitesimal a leaves the 
Hamiltonian H invariant. The only change comes from 
the source term: Hr — > Hy + aJr- This implies: 



= 6Z[T] = -ia / dt (J r (t))r 



(14) 



where (--^r denotes the expectation value in equilib- 
rium with a non-zero T. Since an equilibrium expecta- 
tion value should be time-independent, it follows that 
(Jr(i))r = 0. 

When T -> T(t) the Hamiltonian H ->■ H[T(t)}, and 
the system is driven out of equilibrium. Then the in- 
stantaneous backscattering current can be non-zero. In 
order to calculate this current we introduce ptiie-|generat- 
ing functional for a closed time path (CTP)ll2lll3: 

Z[r+,T_] = (0|S r _(-co,co)S r+ (co,-co)|0), (15) 

where |0) is the state at time t = — oo when the source 
T(t) is adiabatically switched on. The <S-matrix 



Sr^oo, — oo) = Texp 



; f dt [h + H r+ {t)\ 



evolves the state |0) forward in time with a source T + (i), 
while Sr_ (— oo, oo) evolves this state backward in time 



with a source T_(t). Note that choosing T + = T_ for all 
times makes Z[T] = 1 . We also note here that our choice 
of a single initial state is valid only for the zero tem- 
perature formalism. In general a non-zero temperature 
(or any other initial mixture of states) can be straight- 
forwardly accommodated by including an additional S- 
matrix operator S(— oo — i/3, — oo) and tracing over all 
initial statestj. 



1. Ward-Takahashi Identities 

Because of the global chiral symmetry of H, not all 
Green functions generated by the functional -Z[r±] are 
independent. General identities between Green func- 
tions of different orders of the local operator Q can be 
derived by imposing the constraint that the generating 
functional, and any expectation value derived from it, 
is invariant under the global chiral transformation (|^). 
This implies that: 



0=—(Q(t) 

10OL 



Q(t) 



dt'([Q(t),Mt'))), (16) 



where the angular brackets denote expectation value 
of the Heisenbcrg operator at time t, and Jt(t) = 
«r(r)Q(r) + H.c Note that there is no restriction on 
the time dependence of T(t), so that these relations 
constitute non-equilibrium generalizations of the Ward- 
Takahashi (WT) identities derived in equilibrium field 
theoryt3. In terms of operators defined earlier, the WT 
identities are: 



H r (t) 



dt'(jJv(t),Jr(t')]) (17a) 
(17b) 



J T (t)) = H dt'([H r (t),Jr(t')} 

oo 



2. Pumped charge as a consequence of phase change of F 

In order to calculate the pumped charge let us consider 
a time interval [to ,ti], with to < ti, in which the pumping 
parameter can be written as T(t) = To(t) + ST(t), where 
6T(t) vanishes smoothly outside the interval. We can now 
write our Hamiltonian as 7io + STtr , where 5H.r vanishes 
with the perturbation ST(t) outside the chosen time in- 
terval. Then the current generated by backscattering can 
be written in the interaction picture with respect to the 
Hamiltonian H.q: 



I(t) = (S^t,t )Jr(t)S(t,t Q )), 



(18) 



where all operators are in the interaction picture with 
respect-to H , and J r (f) = iT(t)Q(t) + H.c. The S- 
matrixll3 in the interaction picture is: 



S(t,t )=Te- iJ * t ° dt ' SHr{t,) . 



(19) 
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Here T stands for time ordering, and 8Hp{t) is the time- 
dependent perturbation in the interaction representation. 

As a consequence of global chiral invariance, when the 
pumping parameter T(t) has a constant phase <f>o in the 
time interval [ti,tf] (where the system is in an equilib- 
rium state for times t < U) then the pumped charge in 
that interval 



l) 



z [r (+),r(-)]= f tf dt(jr (t (t)) = o, 



(20) 



Therefore, a phase change is the single important pump- 
ing parameter. Notice that in order to achieve a time- 
dependent phase we need to vary a minimum of two ex- 
perimental parameters. 



C. Charge pump: backscattering conductance 

Consider now a change in the phase of T(t) as a per- 
turbation: 



r(t) = |r(t)| e ^ (t) , 

8<j)(t) = 4>{t) - <j)(t ), and 

r (t) = |r(i)|e^ (to) . 



(21a) 



Then, the perturbation parameter can be explicitly writ- 
ten as: 



5T(t)=T (t) cos 8(j>(t) +i sin 8cj){t) - 1 



(21b) 



Therefore, for a small enough change 8<f>(t) in the interval 
[to, t±], we can evaluate ([Lq) perturbativcly to first order 
in 8<t>(t): 



I(t) = (jr (t))-8<t>(t)(H ro (t)) 
-i f dt'8cb(t')([J ro (t),J ro (t')} 

Jt N 



(22a) 



where all operators are in the Heisenberg representation 
with respect to the unpertu rbed Hamiltonian 7i[r ]. 
Using the WT identities ([]7a])-(|l7t]) in (||) we obtain: 



I(t)~(Jr (i))+i / dt'S^t') 



to 



(22b) 



The pumped charge in a small time interval [io^i] can 
be written as: 



x / dt"([J Vo {t),J To (*")]), 



8Q(ti;t ) =i I dt I dt' 8<p(t') 

J to J to 

ft' 



(23) 



x / dt"([J T MJr (t")]), 



where the first terms contribution vanishes due to Q2C 
From the above expression we see that the role of WT 
identities is to ensure that a non-zero pumped charge 
is a consequence of temporal variation in the phase of 
the backscattering amplitude. We note that in obtain- 
ing Eq. (|23|) we have explicitly shown that tp(t) is the 
pumping parameter. 

The approach introduced above can be used to find 
the current at any other time by using the following 
scheme. We divide the entire path traversed by T into 
a sequence of N sub-intervals, as shown in Fig. ||, start- 
ing from the equilibrium state at time t = to with a 
coupling T(to) = r, and returning to the same value at 
tjv = T. The length of the n-th interval [t n ,t n +i] is 
chosen such that (for t £ [t n ,t n +i]) the phase change 
8<j){t) = 4>(t) — 4>(t n ) is small enough for the linear ap- 
proximation to hold. The current I(t) for t £ [t n , t„+i] is 
evaluated by going to the interaction representation with 
respect to the Hamiltonian W[r„] = H + Hr n - Here flr„ 
has a time-dependent coupling T n (t) that has a constant 
phase in the interval t £ [t n ,t n+ i] . We can therefore 
write the current as: 

i(t) = (sHt,t n )j T (t)s(t,t n )) , 

where (24) 

S(t,t n )=Te- lI ^ dt ' 5Hr(t,) . 

Here 8Hr{t) vanishes for t ^ [t n ,t n+ i]. Within this lin- 
ear response approximation we can write the expression 
for the total charge pumped in a cycle as a sum of con- 
tributions arising from each of the intervals: 



iY 



A' 



/ dt i(t) = J2 5Q{t n 



(25) 



Because the various time intervals differ only in the 
value of the parameters that characterize the initial time 
of each interval, it is sufficient to calculate the charge 
pumped in any one time interval. Contributions from all 
other intervals can be obtained by appropriate relabel- 
ing of these parameters. Thus, for example, the c harg e 
pumped in a time interval [to,ti] is obtained from (22t), 
after exchanging the order of t and t' integrations: 

8Q(h) = 1^ dt' 84>{t') r dt f dt" K { o R \t;t") (26) 

J to J V J — oo 

where the retarded backscattering current-current corre- 
lator for t £ [to,ti] is: 

K^(t;t") = ie(t-t")([J To (t), Jr (t")})- (27) 

The charge pumped in any other interval can now be 
found by substituting T «-> r„, and 8T(t) <-> ST n (t). 
Therefore the total charge pumped in a cycle can be ex- 
pressed as: 



Q 



N 

n=0 J t n -i 



dt' 8<j){t') 



dt 



dt" ICi R \t;t"). 



(28) 
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The function multiplying 8<f>{t') is to be thought of as 
a generalized non-equilibrium conductance Q for the 
backscattcring current. We write the total charge 
pumped in a cycle as: 



jv 



= E / dt ' 7-^w«(0> 



2tt 



(29) 



Qn = 2tt / F n (t-t')dt / dt" ICl R \t;t"), (30) 

Jt' J-oo 

where F n (t — t') has been introduced as a smooth cut-off 
function that vanishes as t — > t n . In general, the quan- 
tity Q n is dependent on the precise form of F n . How- 
ever, when the contribution from the end points is neg- 
ligible compared to the contribution from the rest of the 
time interval t £ (i n ,i„_i), we can expect some interval- 
independent behavior. We therefore introduce a simple 
form for the cut-off function, and write: 



e-^-^dt dt" JC^^t"), (31) 



2tt 



where tu n — > + in the aforementioned limit. Next, we 
note that the retarded correlator JC n R ^ is in general not 
time translation invariant, because of the time-dependent 

coupling r(i). The function J^dt" JCi R) (t;t") has a 
maximum along t' = t. For non-interacting electrons, 
because of the singular nature of fermionic Green func- 
tion, the retarded correlator JCn R (t;t") has a maximum 
near t" « t + 1/Ao, where Ao is the upper cut-off of the 
problem, i.e., A < Ep (see Section |lV| for details). Also 
the width of this maximum St* ~ 1/Ao. As a result, 
for a frequency £1 >C A of variation of T, we can treat 
the time-dependent parameters T(t") = T(t) = T(t'), 
with an error of the order of f2/A <C 1. Keeping in 
mind the other energy scale in the problem (the inverse 
length scale of the barrier, Ai , which defines the validity 
of single impurity approximation), the adiabatic limit is 
O <C min[Ao, A^]. In this limit we can replace the time- 
dependent Hamiltonian 7i[r(t)] by the instantaneous (or 
static) Hamiltonian H[T(t')], and treat the retarded cur- 
rent correlator as time translation invariant. 

In the presence of electron-electron interactions in the 
bulk, the nature of singularity of the fermionic Green 
function is modified. This deformation of the Fermi liq- 
uid picture implies that a new energy scale enters into 
the problem, thereby altering the aforementioned adi- 
abatic criterion. In order to identify the correct adia- 
batic limit, we define dimcnsionlcss integration variables: 
t = At; t" = At", whereby the integrand acquires a mul- 
tiplicative factor of A~ 2 , and rewrite Q n as: 

Q n = 2n / e-^di / dt" A~ 2 JC^ ({ + f; - + t'). 
Jo J-oo A A 

(32) 

Because of interactions in the wire, the composite opera- 
tor acquires anomalous dimensions. We write the scaling 



dimension A for the composite operator Q as A = 1 — a, 
where a is the anomalous dimension, and invoke the scal- 
ing hypothesis (for the equilibrium problem scaling holds 
whenever |T| 2 -C 1 for the range < A < Ao; in the non- 
equilibrium case it is further restricted to £1 < A < Ao). 
Since the actual dimension of JC is fixed, we obtain: 



Kf>{t/A + t';i"/A + t') 

= A 2 (A /A) 2a ICi R \t + At'; t" + At')- 



(33) 



As a consequence, the dimensionless coupling constant 

r — > r = (Ao/A) a r, and its time-dependence takes the 

form: r(f2f/A+f2t'). We can therefore identify an energy 
_ l/o 

scale dJr = A r , dependent only on the bare coupling 

r and the upper cut-off Ao. Choosing il < A <C Ao en- 
sures that, for a > 0, the jxnormalizcd (dimensionless) 
coupling constant T 3> The corresponding energy 
scale tup defines a crossover scale that separates two qual- 
itatively different responses to the time-dependence: (i) 
For il <C u>r, the response is an adiabatic modification 
of the ground state with a renormalized coupling T, and 
(ii) for il 3> Wr, the response is a sudden modification 
of the ground state with T = 0. In other words, the sin- 
gular contribution from the retarded correlator JC n (t) 
acquires a width of the order of l/u>r, as compared to 
1/A p for the non-interacting case. For example, in Sec- 
tion VB we show that, for the Luttinger model with in- 
teraction parameter K = 1/2, the singularity of JC^\t) 
at t - 1/Ao has a width St* - l/(A |r| 2 ). Thus, if 
A/uj r < 1, then we can treat T(t") « T(t) w T(t'). 
Therefore, Q n can be approximated by: 



Qn{t') ^ 2ir 



dt 



dt" Ki R \t-t"), (34) 



where Lu n — > 0, because the intervals [t„_i,t n ] have 

lengths - 1/A > 1/cor] iC n R) {t - t") is the retarded 
backscattcring current correlator with the renormalized 
coupling constant r(i'), independent of the integration 
variables t,t". It is important to note that the retarded 
correlator JC is time-translation invariant in the argu- 
ments t, i". This transformation, from a current correla- 
tor without time translation invariance to a renormalized 
correlator with time translation invariance, is explicitly 
shown in Section [v] for the special cases when the scaling 
dimension of Q is 1, and 1/2. A consequence of this time 
translation invariance is that we can rewrite the time in- 
tegrals to obtain: 



Q n ~ 27r / dr 
Jo 



1 



K n R \r) (35) 



The integral involving can be readily identified as 
the conductance for backscattering current at frequency 
iu n . 



T(t') 



27r £i fl) (o + ^»)-£i R) (o) 
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r(t 1 ) = r 1 




r (to) = r 



FIG. 2: Path of the impurity coupling F in the complex plane 



Here (0 + iu> n ) represents the Fourier transform of 
K\i (t) at frequency uj n analytically continued to imagi- 
nary frequency iuj n (with uj n — > + ). Wc emphasize that 
the conductance depends on the renormalized coupling 
T(t') which is time-dependent. As a consequence we can 
write the total charge pumped in a cycle as an integral 
over the path of this coupling constant: 



Q= / dt'j>(t')g m 

J((,(t') = <f>o 



(37) 



Thus the problem of evaluating Q, for f2 <C cur, reduces 
to calculating the dependence of the backscattering con- 
ductance on the renormalized impurity couplings. In the 
next section we shall link this quantity to the more fa- 
miliar conductance of the fcrmionic current. 



D. Pumped charge and dc conductance 



We now proceed to recast the expression (|37[) in terms 
of the response function for fermionic currents j to a 
vector potential A(t): 



G 



5A(t>) 



(38) 



Let us recall that in the case of linear response to a per- 
turbative A which vanishes outside a region of length L, 
the conductance is given by (see, for example, Ref |l4|) 



GM = — 



L/2 



L/2 



dy / dt (1 



[j(o,o)J(y,i)]>. 

(39) 



In order to find an expression for pumped current in 
the low frequency (f2 < up) limit, we begin by writing 
the fermionic current at any point x in the wire at time 



t £ [t n ,t n +i], in a way similar to (|lj) of Section II B, 
albeit now with the effective action with renormalized 
couplings T: 

I(x, t) = (St( tj tn )j(x, t) S(t, t n )) , (40) 

Expanding this perturbatively to first order in the phase 
change 8(j){t) of the coupling T(t) in the ti me interval 
[t n , t n+ {\, and using the notations of Section II B we get: 



I{x,t) = (j(x,t)) 



i / dt'scf>(t')(g(x,t),j tn (t')]). 



(41) 



We can now invoke Eq. (|11| ) and do an integration by 
parts to write the current as: 

I(x,t) = (j{x,t)) 

/•* • r 1/e , - - (42) 

-i / dt'8ct>(t') / dy(\j(x,t),j(y,t')}), 

Jt„ J~l/e 

where e = 2/L < fl. Since the unperturbed Hamiltonian, 
under which all operators evolve, involves T n (t) which 
has a constant phase in the time interval \t n ,t n+ {\, it 
follows that (j(x,t)) = 0. The pumped charge in the 
same interval can be written as: 

SQ(t n ) = - J^ +1 dt'6<p(t') J™ +1 dt K { n R) (t,t') (43) 

where the retarded current correlator: 

K^\t,t')=i dy9(t-t') (fj(x,t)j(y,lf)]), 

and the subscript n reminds us that this depends on the 
time dependent coupling |r n (t)|. In the slow pumping 
limit the contribution to the integral over t from the end- 
point is negligible (or the integral is dominated by 
the singular contribution near t w t'), expression ( |4^ ) 
can be rewritten as: 

dt' 6j>(t') J dt Al fl) (t,i'()14a) 

ntn+l roc 

dt' 5<j)(t') dt[l~e- et ]K^ R) (t,t').(Ub) 



As discussed in section [IC, when fl/cur < 1, we can 
obtain the effective Hamiltonian that evolves the opera- 
tors j. This effective Hamiltonian has static renormal- 
ized couplings T n (t'), yielding a time-translation invari- 
ant correlator: 

Ki R \t,t')->Ki R \t-t'). 



In this limit the expression ( [l4b| ) above, can be identified 
as the equilibrium dc conductance lim c ^ + Gp (0 + ie): 



1 f°° 

G fn (0 + ie)=-- / dt[l-e- et ] 

2 Jo 



71 ' dyi(g(0-,t),j(y,0)Y 



(45) 



where the subscript n reminds us of the renormalizcd cou- 
pling constant T n . On the other hand, expression (44a) 
can be simplified using the current continuity equation 

d x rj(x',t)=d t p{x',t);x'^0, 

inside the expectation value: 



dt 



-l/e 

1 rV« 



dy ([j(x,t),j(y,t')} 



dy 



dx'([d x rj(x',t),j(y,t')] 



dt 

2 « J-i/ 

l/e px—0~ poo 

dy dx' dtd t ([p{x',t),~j{y,t')\ 

-l/e J — oo Jf V 

- / dy dx'([p(x',t'),j(y,t')} 

l/e ,-2:=0 
-l/e 



(_l , )a a ,j( a /-y) = i-. 



(46) 



Here, in the last step, we have used the equal time anoma- 
lous commutator of fcrmionic density and current, 



[p(x,t'),j(y,t')} 



-d x 8(x - y). 



(47) 



Putting all the terms together we get: 



SQ(t„) = I " +1 dt' 6<}>(t') \±- - G f (0 + ie)] (48) 

The prescription to obtain the total charge pump ed in a 
cycle, Q, is exactly the same as in Section [EI Cj : divide 
the entire path of T(t) into appropriate intervals, and find 
the charge pumped in each interval in the linear response 
approximation. We thus find the charge pumped through 
x = in the limit <C ujy'. 



1 



dr r 
T 



1 - lim27rGp(0 + ie) 

e^O 1 



(49) 



III. TRANSPORT OF CHARGE AND SPIN 
USING PARAMETRIC PUMPING 

We now include spin to analyze both the spin pump 
and the charge pump of Ref. |7| We begin by considering 
a Hamiltonian H which conserves the z component of 
the total spin, and chirality in each of the two spin states 
(|, J.). Thus the axial charges 



dx~f' s (x,t)= / dx[jHx 1 t)±jt(x,t)] 



(50) 



V f (t) 



Wire 




FIG 3: Scattering potentials for up (f) and down (j) spins 
whose variation in time drives the spin pump 



%(x,0) 



y (x,o) 




FIG 4: Scattering potentials for up (f) and down (J.) spins 
whose variation in time gives rise to opposite sense of rotation 
for r± 



pump operates by adding to if a time-dependent po- 
tential, locally breaking the chiral symmetry in the spin 
(charge) sector. For the spin pump we consider a local- 
ized magnetic field in the z direction, acting as an impu- 
rity potential V (x,t) that couples anti-symmetrically to 
up (|) and down (j) spin states. In addition we consider a 
localized scattering potential V e {x,t), coupling symmet- 
rically to the spin states. This is shown schematically in 
Figure [| where we introduce the notation = V e ±V . 
The additional term in the Hamiltonian is: 



dx V e (x) tp}- (x)ip<\ (x) + -0 j (x)ipi (x) 
dx V (x) ijjl(x)^(x) - ip\(x)ipi(x) 



(51) 



are conserved. Here j c ' s are the fermion current oper- 
ators for charge and spin currents. The spin (charge) 



As we pointed out before in Section II A, it is the back- 
ward scattering of the right and left movers from the 
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potentials that leads to a non-zero current. Follow- 
ing Ref. we consider the case when the potentials 
V 0<e (x,t) = V ^ e (x)f 0te (t) with real periodic functions 
fo,e(t). In particular, for harmonic pumping we can 
choose f e (t) = cos(f2i+0) and f (i) = cosilt. The elastic 
backward scattering of spin up (down) electrons is then 
given by the complex quantities: 



(52) 



r±(t)= J dxe tM »* [V e (x)f e (t)±V (x)f (t)] 
= V e (2k F )f e (t)±V (2k F )f (t) 



which are the 2fci?-Fourier- transform of the potentials 
that scatter up or down (±) spin particles. Thus, if 
V e (x) and V a (x) are different functions of position (both 
being assymetric with respect to x = 0), then the com- 
plex qu antit ies V e>0 {2kp) have different arguments (cf. 
Section II A ) , and by suitably varying the real functions 
fe,o(t) (making them change sign) we can trace any path 
in the complex plane which winds around the origin. Fur- 
thermore, the the difference between r + and T_ is that 
they wind in opposite directions on the complex plane, 
as shown in Fig. [|. 

We can denote the composite operators, which have 
couplings r±, by Q±. It is more convenient, however, to 
define operators which transform simply under rotations 
in spin space. We therefore write the backscattering term 
in the Hamiltonian (51) by defining linear combinations 

of r ±: 



(53) 



0=o,e 



where the S D ( e ) field is odd (even) under the unitary 
operator lZ y (Tr) = e _? ^ <Ty , which does a 7r-rotation about 
the y-axis in the spin space, and Y e o = ^[r + ± T_]. The 
backscattering current that arises from non-conservation 
of the axial charges J dx j c ' s (x) is now 



J, 



S (c) 



1 



dx d t j s{c) (x) 



(54) 



H.c. 



(o) T «■ e^o(e) 

WT identities, similar to Eq. ( 17a , 17b| ) can now be de- 
rived for each of Q o e fields and their corre spon ding 
back scattering currents. As before (Sections [I C and 
IIP ), WT identities ensure that there is no current with- 
out a change in the phase of r± . 



A. Spin pump 

We now consider the currents generated by a small 
change in the phase of the couplings T± = T e ± r o , in a 
time interval t £ [to,ti\: 



/•< c W)=-(? (c) (x,t))+x; / dtsMt') 



<T=± 



Xi(g s ^(x,t),J! o (t')+aJ? o (t')}) 



(55) 



where, as in Section To(t) is the time dependent 

renormalized coupling with a constant phase for t G 
[t0)*i]- Next, we invoke Eq. (p4]), albeit with renormal- 
ized couplings, and do an integration by parts to write 
the current as: 

i^(x,t) = ]T fdt'\s4> a {t') f dt" 

a= ±Jt 1 J-oo 

x dyi(^ c \x,t),j c (y,t'') + aj s (y,t'' 

J-l/e X 

(56) 

where e = 2/L < f2. The pumped spin (charge) in this 
interval can now be written as: 



5Q s{c) (to)= / dt'J^SMt') 



(57) 



where the retarded current correlator: 



K ( oZM = i I x dy6{t-H) (\j a (x,t),j»(y,t')]y 



Following Section [ID, we rewrite this expression as: 



5Q s{c) (t Q ) = / dt'J2^(t') 



x / e~ et dt 



-\K^ R ], (t)+aK< R ], (t) 

2 |_ 0;s(c),c^ / 0;s(c),sV J 



(58) 



which can be recast as: 
SQ sic) (t )= /V^i^O 

Jto a=± 



x / dt 
lo J 



- \k {R) (t) 

f 1 dt'Y^buf) 

Jtn L Z 



a K 



(R) 

0;s(c),s 



(*) 



(59) 



G s S c) ' c (0 + ie) + aGf c) ' s (0 + ie) 



where 

Gf(0 + ie) 



o 



dt [1 



1/N 



dyi(\j a (0-,t),j b (y,0)] 



(60) 



The mixed conductance G~' c = 0, from the following 
symmetry argument: The Hamiltonian (with the impu- 
rity) is invariant under a combination of spin rotation 
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TLyiji) and local magnetic field inversion, whereas the 
conductance G~' c — * — G~' c . We can now use a less loaded 
notation by labeling the spin (charge) conductance by a 
single superscript. Thus we obtain the following expres- 
sion for total spin and charge transported through any 
point in the wire in a pumping cycle: 



(7— ± 



cr=± 



(61) 



Here G^ s ' ) is the DC charge (spin) conductance with 

renormalized impurity coupling constants T±. For the 
spin pump the choice of time dependence for the param- 
eters V 0j e(x,t) = V ,e{x)f ,e{t) (discussed earlier) implies 
that 



2iri. 



(62) 



We can now write the pumped spin (in integer values of 
h/2) as: 



<T=± 



2m T T 



(63) 
(64) 



B. Charge pump 

The charge pump operates by replacing the magnetic 
field (V (x,t)) by a charge (spin independent) scattering 
potential. The above analysis, and especially the expres- 
sions (^), can be straightforwardly adapted to this case 
by choosing T+(i) = r_(t). We get the charge pumped 
(in units of e) in a cycle: 



ZTTl 

while no spin is pumped. 



-27rG~,, 



(65) 



IV. DISCUSSION: PROBING THE FIXED 
POINTS OF THE ONE DIMENSIONAL 
IMPURITY PROBLEM 

In this section wc discuss the behavior of pumped 
charge and spin for both the spin and charge pumps 
discussed in earlier sections. We find that bulk interac- 
tions play a crucial role in determining what the pumped 
charge or spin will be in both the asymptotic limits of low 
(Q <C u>r) and high (Q 3> u>r) frequency temporal varia- 
tions, where the crossover energy scale up has been intro- 
duced earlier. By contrast, for non-interacting electrons 



there is no such distinction. As we discuss below, it is the 
non- Fermi liquid behavior peculiar to one-dimensional in- 
teracting fermion systems, which is responsible for this 
difference. As a result we find that the pumped charge or 
spin behavior at low pumping frequency probes the fixed 
point behavior of the equilibrium impurity problem. 

The low frequency behavior of pumped charge and spin 
can be obtained from the expressions (37), (|49|), ( p3| ) 
and(|65|), derived earlier. We first consider the case of 
non-interacting electrons. 



A. Non-interacting electrons 



In Section II A we had pointed out that the impurity 
coupling r has zero dimension. In the absence of inter- 
actions this is also the scaling dimension of the renor- 
malized coupling r, which is thereby marginal in the RG 
sense. This absence of scali ng w ith a change of cut-off 
(in the notation of Section II C , the exponent a = 0) 
means that if the barrier V(x) <^ Ep, one can evaluate 
the conductance (of both the fermionic current and the 
backscattering current) at any energy scale using pertur- 
bation theory. To lowest order in the coupling T (a more 
complete calculation is given in Section V A ) we proceed 
as follows. 



1. Free electron Propagator and the Current Correlator 

In order to calculate the retarded backscattering 
current-current correlator )C^[t,t") pcrturbatively in 
the time-dependent coupling T, we need the (local in 
space) free electron propagator for each of the left (right) 
moving fermions. With our normalization (a constant 
density of states per unit length vq = and a cut-off 
A = Ep) this is: 



D(t) = 



it + sgn(i)/A 
The zeroth order term of K,( R > (i, t") is then given by 

Mt-t") 



(66) 



K^(tX) = \T(tW(t")\- 



An 



(t - t"Y + l/Ag. 



(67) 

If the frequency of variation of |r(t")| is much smaller 
than 1/Ao, then we can approximate |r(£")| « |r(i)| in- 
side the integral ft ^dt" t"), as it is dominated 
by the contribution near t' ks t — 1/Aq. Furthermore, 
we can write as |r(i')|. These approximations can 

be performed at every order in the expansion of KS R \ 

so that wc obtain, from J~ dt ft 00 dt"K.( R \t,t"), a DC 
backscattering conductance which depends only the in- 
stantaneous parameter |r(t')|. For the spinless case we 
find the lowest order contribution to the DC backscatter- 
ing conductance: 

Gr^-lrl 2 , 
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Using this expression in Eq. 
charge 

1 



we get a non-quantized 



Q~—(b [T*dT-TdT*]. 

ITT J 

The integral over the path can be rewritten (using Stoke's 
theorem) as an integral over the area of the contour: 



Q 



dA, 



where the differential area element dA is shown in Fig. 
The charge pumped is thus dependent only on the area 
covered in a cycle on the T-plane, not on where the area 
is located in the T-plane. 



B. Interacting electrons: Luttinger liquid behavior, 
spinless case 

In the presence of interactions in the one-dimensional 
system the renormalized coupling T acquires anomalous 
dimension a, which is indicative of the absence of a quasi- 
particle pole in the (bulk) electronic Green function: 

D(k) ~ (|fc| - k F )- 1+a ; 

a = for non-interacting electrons where we recover 
the quasi-particle pole. As a result, the coupling T(A) 
changes with the cut-off A. In the language of the RG, 
the coupling is relevant when a > and irrelevant when 
a < 0: 



dT 
din A 



= aT 



It is wcll-knownlliHi3 that the backscattering impurity 
coupling is always relevant when the electron-electron in- 
teractions are repulsive so that T(A) increases as A — ► fi, 
the frequency of pumping. As £1 — > the DC conduc- 
tance at zero temperature vanishes. This implies that the 
charge pumped in a cycle (^) in the asymptotic limit of 
slow pumping is: 



Q 



l 



dT r 
T 



1 - lim27rGp(0 + ie) 
2iti I T 



(68) 



1. Frequency corrections 

Let us now turn to the question of low-frequency cor- 
rections to the quantization of Q. As £1 approaches lot, 
our assumption of replacing the upper limit of integra- 
tion over variable t in expression (pl[), or (|f3|), by a cut- 
off to n — > + , is no longer accurate. Deviations from this 
arise as the interval length \t n — t n _i\ shortens with in- 
crease of pumping frequency £1 (in order for the linear 



approximation to hold). As a result, the dependence of 
CJp on the end-point of the integration region, i.e., t n , in 
each of the summed quantities in (|8|) grows; likewise for 
G r in ([Is]) . This dependence can be taken into account 
very simply by choosing Lo n = e ~ f2. Then the resulting 
quantity is the conductance G r (0+ie), simply related, by 
analytic continuation, to the real frequency conductance 
G r (e) for e > 0. The functional form G r (e), as e — ► 
is governed by the effective impurity coupling T(A) at 
the energy scale A ~ Q. To determine this behavior we 
have to know the finite ac-bias conductance of the inter- 
acting system with an impurity which requires recourse 
to a particular model. Exact solutions for the Luttinger 
model with a time-independent impurity coupling T have 
been shown to have a scaling form for the conductance 
G(oj) = G{[^-LTb\) with the impurity strength dependent 
energy scale-3 Tb ~ Ao|r| 1 ^ 1_X ' ) , and the Luttinger pa- 
rameter K < 1 for repulsive interactions. In the limit 
wCTb the conductance vanishes as {lo/Tb) 2 ^ k ~ 2 ■ This 
can be used to find the asymptotic frequency dependence 
of the pumped charge as follows. We replace the impurity 
energy scale Tg with the instantaneous one wr(t'), and 
the energy scale set by the applied voltage/frequency, lo, 
by that set by the cut-off e. Thus we obtain: 



Gr(t>)(. € ) 



^r(t') 

The pumped charge is therefore: 
1 f dT 



2/K-2 



(69) 



x Ci£l 2/K ~ 2 \T\ 2/K + C 2 n A ' K ^\T\ A ' K + . 
= 1 - qi (n/u T ) 2 l K - 2 - q^/LOrf'^ + 



(70) 



where C\, Ci are cut-off dependent non-universal con- 
stants, whose ratios are universal numbers. However, in 
the expansion of pumped charge Q the ratios of coef- 
ficients, for example 91/(72, as well as the quantity lot, 
are dependent on the details of the path followed by 
the coupling F and are consequently non-universal. We 
note though, that the pumped quantity (in this case, the 
charge Q) is bounded from above (below) by its value 
in the case of a constant |T| = max(min)[|r(£)|], which 
also corresponds to a larger (smaller) value for the energy 
scale lot ■ The above considerations give us the frequency- 
dependent correction to the quantized pumped charge Q 
in a one-dimensional system of spinless (spin-polarized) 
interacting electrons. 

When the electron-electron interactions are attrac- 
tive, then_.it is known that the impurity coupling is 
irrelevant. As a result the renormalized impurity cou- 
pling at the energy scale A ~ Q, is smaller than the bare 
coupling at energy scale Ao > A and the DC conduc- 
tance G r can be calculated perturbatively. The lowest 
order correction corresponds to setting F = which gives 
zero pumped charge. For the frequency dependent cor- 
rections we turn to the Luttinger model with K > 1. 
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Using the exact resultsB and the scaling arguments used 
earlier for the K < 1 case, we obtain the conductance 
2irGf ~ 1 — q[(n/oj-p) 2K ~ 2 , and the charge pumped in a 
cycle vanishes with a power law: Q ~ (£!/d>r) 2A ~ 2 . 



C. Quantization and fractional charge 

While the above analysis for spinless electrons does 
not directly apply to the case of a quantum wire - where 
both spin species of electrons are present - it correctly 
describes a charge pump operating in a quantum Hall 
system, where the left and right movers are the chiral 
edge excitations. For integer Hall systems the chiral ex- 
citations carry integer charge and correspond to the non- 
interacting case discussed earlier. On the other hand, for 
fractional quantum Hall systems whose edge .excitations 
are described by Wen's edge-state theoryllZHIa, the chiral 
excitations carry fractional charge ve and correspond (in 
our description of the pumping behavior) to interacting 
electrons with K = v = 1/m < 1, where m is the odd- 
integer denominator which characterizes the bulk frac- 
tional quantum Hall state. One interesting question that 
emerges out of the quantization of adiabatically charge in 
units of the electron charge e is how charge fractionaliza- 
tion in strongly correlated one-dimensional systems (in- 
cluding edge state tunneling in fractional quantum Hall - 
FQH - systems) is manifest through pumping. Consider, 
for example, a time-dependent backscattering potential 
in a constriction of a v = 1/3 FQH bar: it will lead to the 
adiabatic transport of an electron charge e within a (slow) 
period t = 27r/f2, or a pumping current I p = ^S7. To 
attain the same current with a non-equilibrium voltage, 

2 

one must have I p = v\V instead. So while the charge 
pumped in a cycle is e, the relationship between the 
pumping frequency S7 and the non-equilibrium voltage V 
is SI = veV/Ti. Both the finite bias and pumped currents 
derive from the phase of a time-dependent backscattering 
potential, and the phase changing rates in the two cases 
are related by the Josephson relation S7 = e*V/h, with 
e* = ve. So even though the charged pumped in a cycle is 
integer, independent measurements of the pumping fre- 
quency and voltage leading to the same current yields the 
fractional charge relation between S7 and V. Simply put, 
the existence of fractional charge in the FQH is directly 
related to the quantization of the Hall conductanccO, 
and to the quantization of pumped charge across the 
edges. 



above in the context of spinless fermions. Assuming that 
interactions in the one-dimensional system respect SU (2) 
symmetry in spin, it is known that the impurity backscat- 
tering coupling r 0ie is relevant whenever the interactions 
are repulsive. As a consequence both the spin and charge 
conductances vanish at low energicstj with power law 
corrections which can be determined by taking recourse 
to a particular interaction model with spin. In particular 
for the Luttinger model one needs to only modify the ex- 
ponents of the renormalization factors in the expressions 
for Q by replacing 1/K — > 1 /K^-llKs and putting 
K s = 2 for the spin isotropic pointtJ ! t§E2l in order to get 
the behavior of Q s ' c . Note that the non-interacting case 
is given by K c = 2, K s = 2, and repulsive interactions 
imply K c < 2. For the spin pump in a wire described by 
the Luttinger model with K s + K c < 2, we obtain: 

2 «(ifc+Trr-i) 



Q' 



l 



n=l 



n—1 v 



si 

Cur ( 



+ 



(71) 



(72) 



For the charge pump the pumped charge in a cycle is 
quantized just as the spin is quantized for the spin pump, 
albeit with different non-universal constants {q^}, as in 
the case of spinless fermions. 

The more general case of pumping in a wire described 
by the Luttinger model with K s ^ 2 can also be stud- 
ied using the general relations derived in Equations ( |63| ) 
and(|65|). All we need, to determine the pumped charge 
and spin in a cycle, is the behavior of the DC conduc- 
tance with renormalizcd impurity couplings. When the 
charge (spin) DC conductance vanishes we have quan- 
tized charge (spin) pumped in a cycle. Allowing for more 
general values of K c and K s we can refer to the plots of 
Kane and Fisher Ref. jl4|] for depicting the regions, in 
the K c — K s plane, where the 2Uf backscattering is rel- 
evant and where it is irrelevant. In the former regions 
the spin pump will transport a quantized spin per cycle, 
and the charge pump will transport a quantized charge 
per cycle, whereas in the latter both will give vanishing 
pumped charge. In passing we would like to note that the 
analysis of parametric pumping presented thus far does 
not apply to the case when more than one backscattering 
operator is relevant. In the Luttinger model this happens 
when K c < 2/9, K s = 2, when both the 2kp backscat- 
tering and 6kp backscattering terms are relevant. 



D. Electrons with spin: spin and charge pump. 

The quantization of the pumped spin in a cycle, Q s , 
depends on the low energy behavior of the DC spin and 



charge conductances G- 



s(c) 



with renormalizcd impurity 

couplings r o e . The behavior of these DC conductances 
can be obtained in a manner similar to that discussed 



E. Finite size and temperature effects 

Real experimental scenarios involve finite length sys- 
tems contacted by wide leads which are usually described 
by Fermi liquid theory. This brings the length of the 
quantum wire Lyy as another important scale in the prob- 
lem. The issue of contacts in the presence of an impu- 
rity in a Luttinger liquid model has been dealt with ex- 
tensively in the literaturdij. Here we use the relevant 
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results. Whenever fl ^S> vp/L\y the (backscattering) 
conductance Qf is determined by the properties of the 
wire alone, so that the above considerations for quantiza- 
tion and corrections from non-zero Q hold. On the other 
hand, when Lw <C vp/Q, Qf is determined by the prop- 
erties of the external leads. In the RG analysis discussed 
earlier, this corresponds to the infrared divergence of T 
being cut off by the energy scale Tivp/ Lyy- Consequently 
the charge (spin) Q c ^ is independent of the pumping 
frequency, and also of the length L w -C vp/VL just as for 
the non-interacting case. To see the finite size scaling one 
can imagine operating the pump in a closed geometry, for 
example, in a ring. Then the finite size effects of Gf. are 
similar to those for the ac-conductance in Ref. |l4| 

The dependence of Q on the DC conductance also im- 
plies that the finite temperature effects (when T > h£l) 
arc similar to that of Ref. [l4| Consequently the low tem- 
perature (fi < T < d>r) expansion is similar to ( |7l|) - 
when T is substituted for f2. 

In conclusion, let ns note here that unlike the case 
of an applied voltageti the frequency corrections to the 
pumped charge are non universal, so that it is not pos- 
sible to write a one-parameter (for example, Tg) scaling 
formula for the pumped charge as a function of frequency. 
This dependence on the shape of a pumping cycle also 
applies for the temperature and finite size effects. 



V. TIME-DEPENDENT IMPURITY IN THE 
LUTTINGER MODEL AT K = 1, ~ 

In this section wc consider electrons with short ranged 
density-density interactions described by the Hamilto- 
nian of a Luttingcr model, written here for the case of 
spinless electrons: 



+ 92 [pR (x)pR (x) + PL (x)p L (x)] | , 



(73) 



with the interactions parameterized by K = 
y(l — 92)1 (1 + .92)- To this wc add the backscat- 
tering Hamiltonian iJr, introduced earlier. 



A. Exact solution for pumped current at K = 1 

In the absence of interactions, 92 = 0, the Luttinger 
parameter K = 1. This case can be readily solved for the 
current. We begin by writing the Hamiltonian with the 
time-dependent backscattering amplitude T(t): 



n 



H + Ht 



= Jdx¥(-ia 3 d x )^(x,t) 

+ s( x ) \r{t)^y R (x,t) + r*(t)Tpy L (x,t)] (74) 



We remind the reader that we have set % = 1 = vp. The 
density of states per unit length for a single spin species 
is isq = 1/(ttKvf) = I/7!", and the constraint that the 
ground state contain a fixed number of particles relates 
the upper cut-off of the theory to the Fermi energy Ep. 
The uniform density of particles in the ground state of the 
wire is kp/ir. In the absence of any time-dependence of 
r the wire is assumed to be in equilibrium. The current, 
which is given by the difference in density of the left 
and right movers, is therefore zero. We now ask for the 
density difference as a function of time at a point in the 
wire to the right of the impurity. For convenience, wc 
choose this point to be in the immediate vicinity of the 
impurity. The time evolution of the fields is (obtained by 
using the Fcrmion anticommutation relations): 

-idt4>R{x) = [H^Rix)] 

= idMx)~6(x)T*(t)i' L {x), (75a) 

-id t tp L (x) = [H,ip L (x)] 

= -tdMx) ~ 5(x)T{t)^ R {x),(75b) 

Since we are interested in the fields in the immediate 
vicinity of the impurity, we integrate these equations to 
obtain: 



A^ R (t) = -iT*{t)\y L {t), 



(76a) 
(76b) 



where 



A^h = tPr(x = 0+) - i/; R (x = 0"), (77a) 
At/>£ = ip L {x = Q-)-ip L {x = + ), (77b) 
4>r,l = TpR, L {x = + )+i> R , L {x = Q-). (77c) 

The field iPr(x, t) being a right-mover is a free field to the 
left of the impurity, not being influenced by the presence 
of the impurity. Likewise for the field ipp(x,t) to the 
right of the impurity. From the relations ( [76] ) we can 
determine the outgoing fields: 

^fl(O-t), (78a) 

j. -r i± wr 

V-L(0 + ,t). (78b) 

The current for free electrons is given by the current in 
the outgoing channel: 



1 


+ |r(t)| 2 


1 


-|r(t)| 2 


1 


+ |r(t)| 2 




-2iT(t) 


1 


+ \m\ 2 


1 


-|r(t)| 2 


1 


+ |r(t)| 2 



~ V (^ (0 ~' *Wx(0-, «)) + H.c, 



(79) 
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as the incoming channel is taken to have zero net current: 

i&£(o+ t)»flnMo + ,t)) - (ip f R (o-,t)id t Mo~,t)) = o. 



Note that use has been made of Eq. (|75j) away from the 
impurity, to obtain time derivatives in the expression for 
the current. 



Using Eq. (78), and the incoming state normalization 
that 

(v4(o-,i)Vii(o- *)} - (^[(o+,<)VL(o + ,t)) = ^, 

we can write the current (at zero temperature) as: 
1 t(t)T*(t)-t*(t)T(t) 



lit) 



^ (i + iiwr 



(80) 



The pumped charge (at zero temperature) can now be 
written as: 



Q 



dt 



1 t(t)T*(t)-t*(t)T(t) 



1 

2?ri 



dT 



(i + |r(i)l 2 ) 

2|r| 2 
(i + |r| 2 ) 2 



c.c. 



(81) 



(82) 



The reflection probability |T| 2 determines the pumped 
charge generated by backscattering of electrons from the 
barrier. The chief feature of non-interacting electron gas 
is the energy independence of the reflection probabil- 
ity (for low-energies). As a consequence, we find that 
the pumped current is determined by the instantaneous 
backscattering amplitude. 



1. Spin pump 

Including spins in the above analysis, to analyze the 
spin pump, is straightforward; different spin species are 
backscattered by different amplitudes (see Section III). 



Therefore, all we have to do to determine the pumped 
charge and spin is to put indices of different spin species 
on r in the expression J80|). The spin current is then 
given by: 

I s (t) = l r (t) - I L (t), (83) 

1 t ± (t)r* ± (t)-t* ± (t)r ± (t) 



(84) 



(i + |r ± (t)py 

The condition for generating a pure spin current requires 
tuning the amplitudes r± such that 7f = — and is 
given by: 

<M*) 



|r+WI 5 



(i + |r + WI 5 



-<M*) 



|r-(f)l s 



(i + |r-WI 2 r 



where (j)±(t) is the phase of r±(i). This requires some 
fine-tuning of the scattering amplitudes. For example, 
in the case of delta function potentials V e , (x) of Section 
III, adjusting the distance t between the barrier and the 
magnetic field to be such that 2kp£ = ir/2 satisfies this 
condition. 



B. Exact solution for pumped current at K = 1/2 

We now consider the case with non-zero interactions 
(.92 7^ 0) in the Luttinger model for spinless electrons. 
With the bulk Hamiltonian H = H LL at K = 1/2, the 
backscattering Hamiltonian H? of Eq. (|?]) has a scaling 
dimension A = K = 1/2 - same as that of a Fermi field 
operator. The low energy behavior of the theory with 
the impurity is therefore identical to that of a system 
of chiral Fermions (different from the original interacting 
fermions) with a HamiltonianEJOti 2 ]: 



dx^(x) [-id x ]i>(x)} 
[A a V'(0) + A* a^t(0) 



i 

7t 



(85) 



where ip is a chiral Dirac Fermion and o is a Majorana 
Fermion representing the impurity. The impurity po- 
tential A = (l/v / 7ra)r, where Ao = 1/ira < Ep is the 
high energy cut-off of the bulk Hamiltonian H, and we 
recognize the dimcnsionless coupling T from the nota- 
tions used earlier. In this section we use the mapping 
to free fermions to find the exact expression for the non- 
equilibrium current arising due to the pumping param- 
eter r — > r(t) acquiring a time dependence. This ex- 
pression was first used in our earlier publication, Ref. [?]. 
Here we show the details of our calculations for the charge 
transported across the wire in a cycle. In the asymptot- 
ically slow limit of pumping (frequency Q/Ep — * 0) the 
pumped charge is shown to be determined by the dc con- 
ductance of the system. Furthermore, a general scaling 
formula for the pumped charge is also conjectured. 
We begin by defining chiral Majorana fermions 



m{x) = — 71 — ; m{x) = *V2 

{rtjix), Vj(x')} = {V(a:), V-V)} = S(x - x 1 )- 
{Vi{x),V2{x')} = ;{nj(x),a} = ;{a,a} = 1, 



(86) 



and denote the real and imaginary parts of the complex 
scattering matrix X(t) as: 

RcA = Ai; ImA = -A 2 ; |A| = A m . (87) 

Next, we write the equations of motion for these fields: 

-id t r]i(x) = [H,r]i(x)] = id x r)i(x) +i\i5(x)a, (88a) 
-id t r]2(x) = [H, r]2{x)] = id x rj 2 {x) + iX 2 S(x)a, (88b) 
*a t a = -[W,o] = iAiTji(0)+iA3Tj2(0). (88c) 

Now we note that the Majorana fields are right-movers 
and therefore, the fields to the right of the impurity are 
dependent on the fields to its left. The latter are, how- 
ever, independent of the impurity potential. 

To find how the chiral Majorana fermion fields at x = 
+ are related to those at x — Q~ , we integrate across 
the impurity, and using the notation 



1 



te(o+) -»&(()-)], 
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obtain: 

-Am = -A m = l -X 2 a; (gg) 

d t a= Ai[A?7i + m(0~)} + + 772(0")]. 

Wc will now define scaled fields: 

Vi -> ^i~V: ^2 -> Aj 1 ^, 
and eliminate a between the equations in (R9J) to obtain: 



-d t Ar]i = 



\\ (Am +771(0-)) +(1^2) 



= -d t An 2 . 



(90) 



Since the scaled fields satisfy A771 = A772, we choose an 
ansatz for relating the fields across the impurity: 



A771 = Mi * jyi + M2 * m = A772 



(91) 



Here * denotes the convolution: 



M\ * 771 



diiMi(i;fi)77i(0-,ii). 



From Eq. (^0|) we get the following differential equations 
for the M's: 



+ ^\p(t-t 1 ) = 0; fori = 1,2 
These can be solved straightforwardly to reveal: 
A|(*i) ( 



(92) 



x cxp 



-0(t-ti) 



dt'-[\\{t') + \l{t')] 



(93a) 



We thus find the relation between the scaled fields on 
either sides of the impurity: 



(93b) 



k=l 



Reverting to original unsealed fields wc define the ker- 
nels: 



M^t-h) = 2M j (t;t 1 )\J 1 (t 1 ), 



(94) 



and write the relations between Major ana fermions across 
x = 0: 



»7i(0+ «) = »h(0- t) + Ai(i) 



Mi * 771 + M 2 * ??2 



(95a) 



77 2 (0+ , t) = 772(0", t) + A 2 (t) Mi * r/f 1 + M 2 * 7? f 2 (95b) 



From the equations (|95j) above we get the relations in 
terms of the Dirac fermions: 



ip(0 + ,t) =7/7(0-, i)+A*(*){Mi * [^ + V f ] 



iM 2 * 



- V] } ■ 



(96) 



The above equation allows us to calculate the current 
at time t. The current is given by the difference in the 
electron density across the impurity 

I(t) = (V>t( +, t)^( +, t) - ^t (0 - t ty( - t) ) i (97) 

where the expectation value is taken in the Hciscnbcrg 
state at a time in the remote past when the system is in 
equilibrium. Using: 



<V^(o-,i)^(o-,t')> = (^(0- t)i&t(o- f)) 

= D(t-t') =^n{E)e iE{ - t - t '\ 

E 

we find: 

I(t) = / df </>(/;') cos[0(i) - 4>{t')} 

J — OO 

x r dt"K^(t,t")- 

j — 00 

where the retarded current correlator 



(98) 



(99) 



cxp 



t" 



nA ° dE 

—n(E) s\nE(t - t") 



in 



2tt 



(100) 



C. Pumped charge: low frequency asymptotics 

Consider the pumping cycle in a time interval [U,tf], 
so that 4>{t) = 0, t < tj. The expression (|99| ) above can 
be integrated over time to find the charge pumped in the 
cycle: 



dt'<j>(t') / d/j cos [0(/j) - </j(t')] 



x / dt"K^{t;t") (101) 
dt'j>(t')g(t'), 



— OO 



We now seek the asymptotic behavior of Q as the ra- 
tio fi/Ao — > 0, where fi is the pumping frequency and 
Aq -Bp is the upper cut-off of the field theory (|85|). 
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Consider first the case of zero temperature. Then the 
energy integral can be performed easily and we obtain: 



K^(t,t") = -\ m ( t )\ m (t")eilf dix ™® 
l-cosA (t-t") 
27r(t - t") 



(102) 



This function has a maximum at t — t" 1/Ao, with a 
width St* « 1/cJr + 1/A ~ l/cj r , where w r = A^(i') < 
Aq. As a result, the time i' has to be within 1/A of 
this maximum, and we can write: X m (Qt) sa A m (fit'). 
We note that the neglected term, 5X m /X m , contributes 
only to even orders to the current, i.e., (SX m / X m ) 2n ~ 

(f2/wr) 2rl , for n = 1, 2 It can therefore be accounted 

for by a suitable renormalization of the energy scale u>r- 
We therefore focus on the analytics of the lowest order 
corrections, and write: 



0(f) 



dt 



A 



dE 
2^ ? 



dt'xKt'yr^ 2 ^'^-^ 



.(E) sm[E(t-t")}. 



(103) 



Using the time-translation invariance of the resulting in- 
tegrand, we have: 



5(f) 



/Aq 7 771 pO 
—n(E) / dr TCj r e^ rT sin[Ei 



L ° dEdn(E) w$(t) 



- - 2G(Q 



_ Ao 2tt d£ cj^(t') + 4E 2 
0). 



(104) 



where G(fi -> 



Luttinger liquidE3c3 with an impurity of strength w 2 



the dc conductance of the K = 1/2 
It has a low temperature (T < loy) behavior G 



(T/^X^t')) 2 . This implies that the pumped charge: 



Q=7T I dt'<j)(t') [1 - 2irG(n -» 0)] 

27T J 

at zero temperature, is Q = e. 



(105) 



by a factor f2): 



dtV(t') / cos [(£(*) -<£(*')] 
it' 

H A m (t)A m (t") 1 ft" dtiA^(ti)/f} 

n 



dt" 

- 30 



x / ^-h(E) sin E(t-t") 

J 27T 



(107) 



From this expression, a large- fl/u>r expansion follows 
when the argument of the exponential is a small quantity, 
i.e., when 



Jq c^iA m 2 (ii) 



O 



«1, 



where X m (t) has its times rescaled according to ( 106 ). 
Expanding the exponential we can compare the resulting 
expansion with the low- frequency expansion of (|70|), by 
substituting for K — * 1/K. We thus find the coefficients 
{q n } for all integer powers n of l/f2 in this dual expan- 
sion. It is clear that unless X m is time-independent, the 
ratios of these coefficients are not universal numbers. 

A particular form of time-dependence for the param- 
eters 4> (t) an d X m (t) based upon the simple model of 
Section flA~| is: 



X 2 m (t) = A 2 1 [l + 5 cos2fit], 0<(5 <1 



(108) 



4>(t) 



tan 



y/l + Sq cos ilt — \/l — So sin fit 



\J\ + So cos ilt + y/1 — So sin fit 



51 



1 + #0 cos 2£lt 



{ L09) 
(110) 



For <5o = this gives a constant and so we can com- 
pute the charge pumped for any frequency fl in terms 
of known functions. We find, after straightforward but 
tedious algebra: 



1 

2^ 



dt'cj>(t') [1 - 2nG(n)} , 



(111) 



where G(£l) is the conduct 
K = 1/2 Luttinger modcll 
is given by: 



ee at a dc voltage Q for the 
^3 and at zero temperature 



D. Crossover temperature and scaling behavior 

In order to look at the high frequency pumping limit, 
we consider the time dependence of <f> and A m with a 
principle frequency f2, so that we can write: 



x m (t) = x m (nt), 4>(t) = <t>(Qt). 



(106) 



Then it is straig htfor ward to see that the expression for 
pumped charge (101) can be written so as to make all 
the fl dependence explicit (upon rescaling all the times 



G(n) = i- - 



An 



2tt 2ttH 



tan 



In Fig. U the pumped charge in a cycle Q is plotted ver- 
sus R 
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= ^ for three values of the parameter 
Sq = 0, 0.5, 0.7. It is clear from the plots that the low 
and high frequency asymptotic behavior is the same in 
all three cases. The uppermost curve is for Sq = - a 
constant A m and therefore a circular pumping cycle and 
has the largest pumped charge Q per cycle for any par- 
ticular value of R. This decrease in pumped charge for a 
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Q 



1 




50 100 150 200 25oR 

FIG. 5: Variation of pumped charge with frequency for three 
different shapes of the pumping cycle 

given frequency is to be expected, as an oscillatory ampli- 
tude leads the response function, Q, out of phase with the 
driving force (j), thereby decreasing the current. As R in- 
creases from 0, deviations from So — curve are apparent 
near R ~ 10; at low frequency (large R) the curves for in- 
creasing S (So — 0.5 and 5 = 0.7 in the figure) converge 
with the curve for So = at progressively greater values 
of R, indicating the decoherent (non-adiabatic) nature of 
deviations from the circular shape of the pumping cycle. 

Thus, although there is a range of frequency where 
pumped charge is dependent on all the details of the time- 
dependent parameters in the problem, the asymptotic 
behavior is universal. A scaling formula for the pumped 
charge can therefore be written when T(t) = T(p.t): 

Q = Q(n/u>r,Scj r /u>r) (112) 

where ujt = § J Q 2 ^^ Q dtX^t) is a cross-over energy scale 
in the time-dependent problem, and Su>r represents a col- 



lection of other (time-independent) p aram eters contained 
in r similar to the quantity So in Eq. ( | | ). The behavior 
of the scaling function is described by: 

Q(0,x) = l, Q(x,y)<Q(x,0). 

Let us note that it has not been proven that the inequal- 
ity will hold for a more general form for X(t). 

VI. CONCLUSION 

In conclusion, we have studied the adiabatic limit of 
a quantum pump in a quantum wire, by mapping it to 
a problem of a time-dependent backscatterer in a Lut- 
tingcr liquid. We have shown that the properties of scale 
invariance and chiral symmetry of the Luttinger liquid 
imply a relationship between the pumped quantities and 
the dc conductance. The difference between a Fermi liq- 
uid and a Luttinger liquid behavior in a one-dimensional 
system is clearly brought out by this mechanism of charge 
and spin transport. It is found that for a Fermi liquid 
the adiabatic regime for which the pumped quantities 
per cycle are independent of pumping frequency - ex- 
tends up to pumping frequency fl < Ep, whereas for 
the Luttinger liquid this regime of pumping is limited to 

<C lut <C Ep. The barrier energy scale wr is found 
to be a crossover energy scale between adiabatic and im- 
pulsive (sudden) response to the time-dependence. This 
distinction is absent in the Fermi liquid picture. An- 
other way of interpreting these results is to associate adi- 
abaticity, for this mechanism of pumping, with the scale 
invariant fixed points of a wire with an impurity. For the 
non-interacting case the beta function for the impurity 
coupling is zero, consequently any pumping frequency 
implies an adiabatic response. For the interacting Lut- 
tinger liquid the low-energy fixed point corresponds to 
an infinite barrier. Therefore, the adiabaticity criterion 
necessarily depends on the ratio fl/uir- 
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Note that this does not mean that the continuum field 
theory itself has a cut-off vf/£, but only that at higher 
energies the potential structure gets resolved. 
Besides the local backscattering operator Q that we explic- 
itly consider, there are other operators present in the bare 
Hamiltonian, that correspond to multi-particle backscat- 
tering processesE-3. In the bare Hamiltonian each of these 
operators, Q' , has a different scaling dimension, A', and a 
different (dimensionless) bare coupling constant, T' , which 
involves 2nfcF-Fourier-transform (n > 1) of the potential 
V(x). As should be clear from the analysis presented here, 
it is only the most relevant operator that determines the 
conductance in the adiabatic limit. 



